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Abstract. The preliminary results of a search for time-integrated CP violation in D° — » h~h + (h = K, n) decays 
performed with 0.6 ftr 1 of data collected by LHCb in 2011 are presented. The flavour of the charm meson is 
determined by the charge of the slow pion in the D* + — > D°n + and D*~ — » D a n~ decay chains. The difference 
in CP asymmetry between D° —> K K + and D° —> n~n + , AA C p = A C p(K~K + ) - A C p(n~n + ), is measured to be 
AA CP = [-0.82 ± 0.21(stat.) ± 0.11(syst.)] %. This differs from the hypothesis of CP conservation by 3.5o\ 



1 Introduction 



The time-integrated CP asymmetry, Acp(f), for the final 
states / = K K + and / = Ti~n + has two contributions: an 
indirect component (to a good approximation universal for 
CP eigenstates in the Standard Model) and a direct com- 
ponent (in general final state dependent). In the limit of U- 
spin symmetry, the direct component is equal and opposite 
in sign for K K + and n ffl- In the St andard Model, CP 
violation is expected to be small [2 3 1]. However, in the 
presence of physics beyond the Standard Model the rate of 
CP violation could be enhanced HI 141 . No prior evidence of 
CP violation in the charm sector has been found. 

The most precise measurements to date of the time- 
integrated CP asymmetries in D° — ♦ K K + and D° — > 
7i -n + were made by the CDF, BABAR, and BELLE col- 
laborations 0516171 . LHCb has previously reported prelimi- 
nary results based on 37 pb~' of data collected in 2010 [8 1. 
In the limit that the efficiency for selected events is in- 
dependent of the decay time, the difference between the 
two time-integrated asymmetries, AAqp = Acp(K~K + ) - 
Acp{n~n + ), is equal to the difference in the direct CP asym- 
metry. However, if the dependence of the efficiency on 
the decay time is different for the K K + and 7r~7r + final 
states, a contribution from indirect CP violation remains. 
The Heavy Flavor Averaging Group (HFAG) has combined 
time-integrated and time-dependent measurements of CP 
asymmetries taking account of the different decay time ac- 
ceptance to obtain world-average values for the indirect 
CP asymmetry of c$ = (-0.03 ± 0.23)% and the differ- 
ence in direct CP asymmetry between the final states of 
Aa% = (-0.42 ± 0.27)% 0. 

In these proceedings, LHCb results 11011 11 for the mea- 
surement of the difference in integrated CP asymmetries 
between D° — > K K + and D° — > 7r~7r + , performed with 
approximately 0.6 ftr 1 of data collected in 2011, are pre- 
sented. The initial state (Z)° or D°) is tagged by requiring 
a D* + — ♦ D°jt + decay. The use of charge-conjugate modes 
is implied throughout, except in the definition of asymme- 
tries. 
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2 Formalism 

The raw asymmetry for tagged D° decays to a final state / 
is given by A law (/), defined as: 

_ N(D* + -> D\f)n + ) - N(D*- -» D Q (/>T) 

Arm/Kj) — ,, — „ - > 

N(D* + -> £>°(/> + ) + N(D*- -> D°(f)7T) 

(1) 

where N(X) refers to the number of reconstructed events 
of decay X after background subtraction. 

The raw asymmetries may be written as a sum of var- 
ious components, coming from both physics and detector 
effects: 

A raw (/) = Acp(f) + A D (f) + A d (tt s ) + A P (D* + ). (2) 

Here, Acp(f) is the intrinsic physics CP asymmetry, Aq(/) 
is the asymmetry for selecting the D° decay into the fi- 
nal state /, Ad(tt s ) is the asymmetry for selecting the 'slow 
pion' from the D* + decay chain, and A P (D* + ) is the produc- 
tion asymmetry for prompt D* + mesons. The asymmetries 
Acp, Ad and Ap are defined in the same fashion as A raw . 

For a two-body decay of a spin-0 particle to a self- 
conjugate final state there can be no D° detection asym- 
metry, i.e. A D (K~K + ) = A D (^7r + ) = 0. Moreover, to first 
order Ad(tt s ) and Ap(D* + ) cancel out in the difference 

A mvj (K-K + ) - A raw (7T7r + ), 

leaving a quantity, defined as AAcp, equal to the difference 
in physics asymmetries: 

AAa- = A CP (K-K + ) - A CP (7T-n + ), (3) 
= A raw (/r/T) - A raw (7T7r + ). (4) 

To minimize second order effects, related to the slightly 
different kinematic properties of the two decay modes, the 
analysis is done in bins of the relevant kinematic variables, 
as shown later in Sees. [4] and [5] The physics asymmetry of 
each final state may be written at first order as lfT2ll : 

Acp(f) * a%(J) + ^fl™ d , (5) 



EPJ Web of Conferences 



where a^pif) is the asymmetry coming from direct CP vi- 
olation for the decay, (t) is the average decay time in the 
sample used, r the true D° lifetime, and a 1 ^ is the asymme- 
try associated with CP violation in the mixing. To a good 
approximation this latter quantity is universal [ 1 1, and so 

AA CP = [a%(K-K + ) - ag(>f> + )] + ^fl'^ d , (6) 

where A{t) is the difference in average decay time of the 
D° mesons in the K K + and n ~n + samples. In the limit 
that A{t) vanishes, AAcp probes the difference in direct CP 
violation between the two decays. 



3 Dataset and selection 

A description of the LHCb detector may be found in Ref. ITPJI . 
The field direction in the LHCb dipole is such that charged 
particles are deflected in the horizontal plane. The current 
direction in the dipole was changed several times during 
data taking; about 60% of the data was taken with one po- 
larity and 40% with the other. 

Selections are applied to provide samples of D* + — > 
D°7T + candidate decays, with D° -> K~K + and 7T7r + . A 
loose D° selection including a mass window of full width 
100 MeV/c 2 was already applied during data taking, in 
the final stage of the high level trigger (HLT). In the of- 
fline analysis only candidates that were accepted by this 
trigger algorithm are considered. Both the offline and HLT 
selections impose a variety of requirements on kinemat- 
ics and decay time to isolate the decays of interest, in- 
cluding requirements on the track fit quality, on the D° 
and D* + vertex fit quality, on the transverse momentum of 
the D° (pr > 2 GeV/c), on the decay time t of the D° 
(ct > 100 |^m), on the helicity angle of the D° decay, that 
the D points back to a primary vertex, and that the D° 
daughter tracks do not. Fiducial requirements are imposed 
to ensure that the tagging soft pion lies within the central 
region of the detector acceptance. In addition, the offline 
analysis exploits the capabilities of the RICH system to 
distinguish between pions and kaons when reconstructing 
the D°. 

Defining the mass difference as 5m = m(h + lTn + ) - 
m(h + h~) - m(n + ), the mass and mass difference spectra 
of selected candidates are shown in Figs. [T] and [2] respec- 
tively. The difference in width between the K K + and n ~n + 
mass spectra arises from the different opening angles of 
the two decays. The D* + signal yields are approximately 
1.44x 10 6 in the tagged K~K + sample, and 0.38x 10 6 in the 
tagged n~n + sample. The fractional difference in average 
decay time of D candidates passing the selection between 
the K~K + and iCi? samples is A(t)/ T = (9.8 ± 0.9)%. 



4 Fit procedure 

Fits are performed on the samples in order to determine 
A lm (K~K + ) and A raw (7r~7r + ). The analysis is performed in 
54 kinematic bins, divided by pj and pseudorapidity (77) 
of the D* + candidates, momentum of the tagging soft pion, 
and whether the initial trajectory of the slow pion is to- 
wards the left or right half of the detector. A binning is 
imposed for the reason that the production and detection 
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Fig. 1. Fits to the m(K~K + ) and m{n~n + ) spectra of D* + candi- 
dates passing the selection and satisfying < dm < 15 MeV/c 2 . 
The vertical blue lines indicate the signal window of 1844- 
1884MeV/c 2 . 



asymmetries can in general vary with p T and 77, and so can 
the detection efficiency of the two different D° decays, in 
particular through the effects of the particle identification 
requirements. 

The events are further partitioned in two ways. First, 
the data are divided between the two dipole magnet set- 
tings. Second, the first 350 pb~' of data are processed sep- 
arately from the remainder, with the division aligned with 
a break in data taking due to a LHC technical stop. In total, 
therefore, 216 independent measurements are considered 
for each decay mode. 

One-dimensional fits to the mass difference spectra are 
performed. The fits include separate components for the 
signal and background lineshapes. The signal is described 
as the sum of two Gaussian functions with a common mean, 
convolved with an asymmetric function. The background 
is described by an empirical function of the form 

h _ e -(ftn-ftno)/cj 

where cw«o and c are parameters describing the threshold 
and shape of the function, respectively. In each case an un- 
binned maximum likelihood fit is used. The D* + and D* 
samples are fitted simultaneously and share several shape 
parameters, though a charge-dependent offset in the cen- 
tral value and overall scale factor in the mass resolution 
are allowed. The raw asymmetry in the yields of the sig- 
nal component is extracted directly from this simultaneous 
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Fig. 2. Fits to the mass difference spectra, where the D° is recon- 
structed in the final states K K + (top) and rTn* (bottom), with a 
D° mass lying in the window of 1844-1884 MeV/c 2 . 

fit. An example fit from one measurement bin is shown in 
Fig- i 

The one-dimensional fits used for the tagged data do 
not distinguish between correctly reconstructed signal and 
backgrounds that peak in the mass difference. Such back- 
grounds can arise from D* + decays in which the correct 
slow pion is found but the D° is partially mis-reconstructed. 
However, these backgrounds are suppressed by the use of 
tight particle identification requirements and a narrow D° 
mass window. From studies of the D° mass sidebands this 
contamination is found to be approximately 1 % of the sig- 
nal yield. 



5 Results and systematic uncertainties 

A value of AAqp is determined in each measurement bin us- 
ing the results for A rdvJ {K~ K + ) and A raw (7r~7r + ). The x 2 /ndf 
of these measurements has a value of 21 1/215. A weighted 
average is performed to yield the result AAqp - (-0.82 + 
0.21)%. 

Numerous robustness checks are made, including mon- 
itoring the value of AAcp as a function of time (Fig.|4ji, re- 
performing the measurement with more restrictive RICH 
particle identification requirements, and using a different 
D* + selection. Potential biases due to the inclusive hard- 
ware trigger selection are investigated with the subsample 
of data in which one of the signal final-state tracks is di- 
rectly responsible for the hardware trigger decision. In all 
cases good stability is observed. 



Fig. 3. Example fit used in the AAcp analysis. The first kinematic 
bin of the first run period with magnet up polarity is shown for 
the D° -> K K+ final state. 
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Fig. 4. Time-dependence of the measurement. The data are di- 
vided into 19 disjoint, contiguous, time-ordered blocks and the 
value of AA C p measured in each block. The red dashed line shows 
the result for the combined sample. 



Table 1. Summary of absolute systematic uncertainties for A A C p. 



Effect Uncertainty 

Fiducial requirement 0.01% 

Peaking background asymmetry 0.04% 

Fit procedure 0.08% 

Multiple candidates 0.06% 

Kinematic binning 0.02% 

Total 0TT% 



Systematic uncertainties are assigned by loosening the 
fiducial requirement on the soft pion; by assessing the ef- 
fect of potential effect peaking background in toy Monte 
Carlo studies; by repeating the analysis with the asymme- 
try extracted through sideband subtraction instead of a fit; 
with all candidates but one (chosen at random) removed in 
events with multiple candidates; and comparing with the 
result obtained with no kinematic binning. In each case the 
full value of the change in result is taken as the systematic 
uncertainty. These uncertainties are listed in Table [T] The 
sum in quadrature is 0. 1 1 %. 
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6 Conclusions 

LHCb has measured the time-integrated difference in CP 
asymmetry between D° — > K K + and D° — > n~n + decays, 



AA 



CP 



= ( 



a-(K-K + ) - a™(n-n + )) + 0.098 a™*, to be 



,dir/ 



AA CP = [-0.82 + 0.21(stat.) ± 0.11(syst.)] % 

with 0.6 fbr 1 of 2011 data, where the first uncertainty is 
statistical and the second systematic. Combining the statis- 
tical and systematic uncertainties in quadrature, the signifi- 
cance of the measured deviation from zero is 3.5<x. The re- 
sult is consistent with the current HFAG world average [9 1. 

Subsequent to this result being presented, several pa- 
pers on the theoretical interpretation have been written, ex- 
amining the possible size and uncertainty in the SM contri- 
bution and considering possible new physics contributions. 
A partial list may be found in Ref . H14I15I16I17I18I19I20I21I221 
We thank our theory colleagues for their help in under- 
standing this result. 
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